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Efficient recognition of totally nonnegative matrix cells 

S Launois* and T H Lenagan 



Abstract 



The space oi m x p totally nonnegative real matrices has a stratification into totally 
nonnegative cells. The largest such cell is the space of totally positive matrices. There is a 
well-known criterion due to Gasca and Pefia for testing a real matrix for total positivity. This 
criterion involves testing mp minors. In contrast, there is no known small set of minors for 
i^H ■ testing for total nonnegativity. In this paper, we show that for each of the totally nonnegative 

■ cells there is a test for membership which only involves mp minors, thus extending the Gasca 

and Pefia result to all totally nonnegative cells. 



o: 

^ ■ 1 Introduction 
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An m X p matrix M with entries from M is said to be totally nonnegative (tnn for short) 
if each of its minors is nonnegative. Further, such a matrix is totally positive if each of 
its minors is strictly positive. (Warning: in some texts, the terms totally positive and 
strictly totally positive are used for our terms totally nonnegative and totally positive, 
respectively.) 

Totally nonnegative matrices arise in many areas of mathematics and there has been 
considerable interest lately in the study of these matrices. For background information 
and historical references, there are the newly published books by Fallat and Johnson, [3] 
and Pinkus, and also two good survey articles [1] and [5]. 

In general, the number of minors of an m x m matrix is XlfcLi (T)^ ~ (nT) ~ 1 ~ '^^^ 
by using Stirling's approximation to the factorial. Thus, it is obviously impractical to check 
whether a matrix is totally nonnegative or totally positive by naively checking all minors. 



*The first author is grateful for the financial support of EPSRC first grant EP/I018549/1. 
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In this paper, we are concerned with the recognition problem: can one decide, by check- 
ing a small number of minors, whether or not a real matrix is totally nonnegative or totally 
positive? There is an efficient set of minors to test for total positivity: Gasca and Pena, 
[6] , specify mp minors that can be used to check total positivity. In order to describe their 
result, we need to introduce some terminology. If a minor of a matrix is formed by using 
rows from a set I and columns from a set J then we denote the minor by or [/| J](M) 

if it is necessary to specify that we are using the matrix M. A minor is said to be an initial 
minor if I and J each consists of a sequence of consecutive integers and 1 G / U J. Each 
initial minor is specified by its bottom right entry and so there are mp initial minors in an 
m X p matrix. Gasca and Pena's criterion is that a real matrix is totally positive if each 
of its initial minors is strictly positive. The number of minors used in Gasca and Pena's 
criterion, mp, is best possible. However, there are many other possible choices for a set of 
mp minors which will give an efficient check for total positivity. Interested readers should 
consult [5] where the notion of double wiring diagrams is discussed and it is shown that 
each double wiring diagram gives rise to a set of minors that can be used to check for 
total positivity of n x n real matrices, see [5i Theorem 16]. 

In contrast, as far as we know, there is no small set of minors that can be used to check 
for total nonnegativity. The best result we are aware of is a result by Gasca and Pena, [7] , 
which specifies a set of minors of size approximately the square root of the total number of 
minors of the matrix that can be used to check total nonnegativity of an invertible matrix, 
see [SI Theorem 29] and the remark following that theorem. 

This discrepancy in the two tests may be explained by considering the notion of the cell 
decomposition of the space of totally nonnegative matrices that was studied by Postnikov 
in [13]. 

Let Ai = Ai^piM) denote the set of m x p matrices that are totally nonnegative. 
Let Z he a set of minors and denote by the set of matrices M G A^^°p(]R) for which 
[/|J](M) = if and only if [I\J] G Z; this set is known as the (totally nonnegative) cell 
corresponding to the set Z. For many choices of Z, the cell corresponding to Z will be 
empty: for example, a 2 x 2 matrix has five minors, so there are 32 choices of Z but it 
is not too difficult to check that only 14 of these choices give rise to nonempty cells. For 
example, it is an easy exercise to check that the totally nonnegative cell corresponding to 
the choice Z := {[1|1]} is empty. Nevertheless, in general, there are still a large number 
of non-empty cells: a weak, but easily established, lower bound is that A^^°p(R) has more 
than m^ nonempty cells. 

The nonempty cells provide a stratification of the set A4. The cell corresponding to 
Z = is the set of totally positive matrices. Although one can show that the set of totally 
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positive matrices is dense in the set of totally nonnegative matrices, nevertheless, it is only 
one of many nonempty cells. 

The main aim of this paper (achieved in Theorem I4.4p is to provide a set of mp minors 
for deciding whether or not an arbitrary real mxp matrix belongs to a specified nonempty 
cell. This result extends the criterion of Gasca and Peha for the totally positive cell to any 
nonempty cell. In the special case of invertible square matrices, our result also extends 
related work of Fomin and Zelevinsky. In [H Theorems 4.1 and 4.13] Fomin and Zelevinsky 
give recognition criteria for those totally nonnegative cells in n x n square matrices which 
do not contain the determinant. Note however that there are many totally nonnegative 
cells in A^^°(M) that do contain the determinant. For example, in the 3x3 case, there 
are 230 nonempty totally nonnegative cells, of which 194 contain the determinant, and in 
the 4x4 case, there are 6902 nonempty totally nonnegative cells, of which 6326 contain 
the determinant. 

The methods used to construct our test were originally developed by Cauchon, [2], 
while studying quantum matrices. In recent papers, a close connection has been shown 
to exist between the cell decomposition of totally nonnegative matrices and the so-called 
invariant prime spectrum of quantum matrices, see, for example, [HIEIIII] and it has become 
apparent that methods used in the quantum world may be used to suggest methods and 
results in the classical world of total nonnegativity, see, for an example, [10]. 

The tools we use here depend on the notion of Cauchon diagrams and the deleting 
derivation algorithm that Cauchon developed in [2] for quantum matrices. Here, we pro- 
pose that the corresponding algorithm for totally nonnegative matrices shall be called the 
Cauchon reduction algorithm. 

In the next section, we describe Cauchon diagrams and the Cauchon reduction algo- 
rithm. We also introduce results from [8] which show how mxp Cauchon diagrams can be 
used to parameterise the nonempty totally nonnegative cells in A^^°p(M). As an example, 
we consider the totally positive cell. This example motivates the discussion of lacunary 
sequences for Cauchon diagrams which are introduced in Section [31 Finally, in Section HJ 
we show how, for each Cauchon diagram, lacunary sequences can be used to specify a set 
of mp minors which can be used to test for membership of the totally nonnegative cell 
corresponding to the given Cauchon diagram. 

2 Cauchon diagrams and Cauchon reduction 

Set |1, m] := {1,2,..., m} for any m > 1. 
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Definition 2.1. |^ An m x p Cauchon diagram C is simply an m x p grid consisting of 
mp squares in which certain squares are coloured black. We require that the collection 
of black squares has the following property: if a square is black, then either every square 
strictly to its left is black or every square strictly above it is black. 
We denote by Cm,p the set oi m x p Cauchon diagrams. 

Note that we will often identify an m x p Cauchon diagram with the set of coordinates 
of its black boxes. Indeed, if C G Cm,p and (i, a) G |1, m] x we will say that (i, a) E C 

if the box in row i and column a of C is black. For instance, for the Cauchon diagram 
C of Figure [H we have (2, 4) G C, whereas (4, 3) ^ C . (Note that we are using the usual 
matrix notation for the (i, a) position in a Cauchon diagram; so that the (1, 1) box of a 
Cauchon diagram is at the top left.) 
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Figure 1: An example of a 4 x 5 Cauchon diagram 



Definition 2.2. Let X = { real m x p matrix and let C be a Cauchon diagram 

(of size m X p). We say that X is a Cauchon matrix associated to the Cauchon diagram C 
provided that for all (^,a) G x we have Xi^a = if and only if {i,a) G C. If 

X is a Cauchon matrix associated to an unnamed Cauchon diagram, we just say that X 
is a Cauchon matrix. 

A key link between Cauchon diagrams and totally nonnegative matrices is provided by 
the following easy lemma. 

Lemma 2.3. Every totally nonnegative matrix overM. is a Cauchon matrix. 

Proof. Let X = ( tnn matrix. Suppose that some Xi^a = 0, and that Xk^a > 

for some k < i. Let 7 < a. We need to prove that Xj^^ = 0. As X is tnn, we have 
—Xk,aXi,'y = det ( Xi,'^ Xi,'" ) > 0. As Xk,a > 0, this forces Xi^^ < 0. But since X is tnn, we 
also have Xi^^ > 0, so that Xi^-y = 0, as desired. 

Therefore X is a Cauchon matrix. □ 

In order to define the Cauchon reduction algorithm we need the following notation. 
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Notation 2.4. • We denote by < the lexicographic ordering on N^. Recall that 

{i, a) < U, 13) [{i < j) or {i = j and a < /?)]. 

• Set E° = (|l,ml X \ {(1, 1)} and E = E° U {{m + l,p)}. 

• Let {j,/3) € E°. Then denotes the smallest element (relative to <) of the set 

{{i,a) G E I < {t,a)}. 

Thus, = {j,/3 + 1) when /3 < p and = (j + l,p). 

Convention 2.5 (Cauchon reduction algorithm). Let M = (xj^^) be a real m x p matrix. 
As r runs over the set E, we define matrices M'^*'-' := fx,-^2) as follows. 



t,aJ 



1. When r = (m + the entries of the matrix M^"^^^'^^ are defined by 



X. 



{m+l,p) 



Xi^a for all {i,a) e [l,m] x so M^'^+^'P) = M. 



2. Assume that r = (j, /3) G -E° and that the matrix M*^*"^) = (Xj-*^^) has already been 
defined. The entries x-'^ of the matrix M*-^^ are defined as follows. 

(a) If xJ.J^ = 0, then xj^] = xj^^^^ for all {i,a) G |l,ml x 

(b) If x'^'P ^ and {i,a) G |l,m] x then 

„(0 _ J ^ - ^1 (4?) 4i ^ if z < J and a < /3 

■* otherwise. 



We say that M^^^ is the matrix obtained from M by applying the Cauchon reduction 
algorithm at step r, and Xj ^ is called the pivot at step r. 

3. If r = (1,2), then we set ti^a '■= x]^ for all {i,a) G x |1,J9]. Observe that 



x^ll = X •''q, ■* for all r < {i,a), and so ti^a = xf'^^ = xf'"^ for all {i,a) G E°. 
The matrix M := M^^''^^ is t/ie matrix obtained from M at the end of the Cauchon 
reduction algorithm. 

In [8] it is proved that the Cauchon reduction algorithm provides a bijection between 
the nonempty cells of Ai^p{M.) and the Cauchon diagrams of size m x p. Let us give more 
details about this bijection. First, let us recall the following result from fS]. 

Theorem 2.6. A real mx p matrix M is totally nonnegative if and only if M, the matrix 
obtained from M at the end of the Cauchon reduction algorithm, is a nonnegative Cauchon 
matrix. 
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Proof. This is the content of [H Theorem 4.1 and Corollary B5]. □ 

Thus, each totally nonnegative matrix is, via Cauchon reduction, associated to a Cau- 
chon diagram. In other words, we have a mapping n : M C from Ai^.p{^) to the 
m X p Cauchon diagrams, where C is the Cauchon diagram associated to the 
Cauchon matrix M. Moreover, if M and N are two tnn matrices such that 7r(M) = vr(A^), 
then it follows from [8t Corollary 3.17] that M and belong to the same tnn cell. So, each 
nonempty tnn cell in A1^*^p(M) is a union of fibres of it. On the other hand, it follows from 
[T3] (see also [H Section 6]) that the number of nonempty tnn cells in A^^°p(]R) is equal to 
the number of m x p Cauchon diagrams. Thus, the nonempty tnn cells are precisely the 
fibres of vr, and we have just proved the following result. 

Theorem 2.7. 1. Let M and N be two tnn matrices. Then M and N belong to the 
same tnn cell if and only if tt{M) = vr(A^), ie if and only if M and N are associated 
to the same Cauchon diagram. 

2. Nonempty tnn cells are parametrised by Cauchon diagrams, and the nonempty tnn 
cells in A^.^°p(R) are precisely the sets 

5° := {M e M^iR) I 7r(M) = C}, 

where C runs through the set of m x p Cauchon diagrams. 

3. Let M G Ai^piM.), and denote by M = (ti^a) the matrix obtained from M at the 
end of the Cauchon reduction algorithm. Then M & if and only if ti^a = z/ 
(i, a) G C and ti^a > z/ {i, a) ^ C . 

Proof. The first two statements follow from the above discussion. The third assertion 
follows from the fact that M is a nonnegative Cauchon matrix associated to C if and only 
if ti^a = if (i, a) e C and ti^a > if {i, a) ^ C. □ 

Let C be a Cauchon diagram. The tnn cell is called the tnn cell associated to the 
Cauchon diagram C. 

As an example, we will consider the totally positive cell. In order to reconcile the 
Gasca and Pena approach with our approach, we need to describe a symmetry of totally 
nonnegative matrices. The reason is that we are using Cauchon reduction, which starts at 
the bottom right of a matrix and moves along columns from right to left, whereas Gasca 
and Pena use techniques such as Neville elimination, which start at the top left and move 
down columns. This makes no essential difference, except to the notation: informally, 
reflection in the antidiagonal preserves total nonnegativity. 
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To be precise, let M = (aij) he an m x p matrix and let iVf be the pxm matrix defined 

by {MP)ij := 

am+i-j,p+i-j- Let / and J be row and columns sets for M^. Then standard 
properties of determinants show that [J| J](M'') = [m + 1 — J\p + 1 — J](M). It follows 
that M is totally nonnegative if and only if is totally nonnegative, and, similarly, M 
is totally positive if and only if is totally positive. (Here, x + I := {x + i \ i E /}.) 

We say that a minor [I\J] is a final minor of an m x p matrix if each of / and J consists 
of a consecutive set of integers and either m E I oi p E J . It is clear that initial minors of 
M correspond to final minors of Af; so the Gasca and Pena criterion can be reformulated 
as: an m X p matrix M is totally positive if and only if each of its mp final minors is strictly 
positive. 

If the Cauchon reduction algorithm is carried out on a totally positive matrix M then 
one can calculate that tj^Q,, the (z, a) entry of M, the matrix obtained at the end of the 
Cauchon reduction algorithm, is given by 

U,a = + 1, . . . , z + r I a, a + 1, . . . , a + r](M) ■ [z + 1, . . . , i + r | a + 1, . . . , a + r](M)~\ 

where r = min{m — i,p — a}. Note that this formula only involves final minors of M. In 
particular, M is a strictly positive matrix and so the Cauchon diagram corresponding to 
the totally positive cell has all boxes coloured white. Also, one can then easily calculate 
that for any final minor + . . . ,i + r \ a, a+, 1 . . . ,a + r ](M) = ti^ati+i,a+i ■ ■ ■ ti+r,a+r- 
In considering totally nonnegative cells other than the totally positive cell, two problems 
arise: some of the ^ will be equal to zero, and some of the (final) minors of the matrix 
M will be equal to zero, and so cannot be inverted. In seeking a criterion for recognition 
of these cells that is similar to the Gasca and Pena criterion for total positivity we were 
led to consider certain sequences called lacunary sequences which were first introduced by 
Cauchon in his study of quantum matrices, |2] . These sequences are described in the next 
section. 

3 Lacunary sequences 

Definition 3.1. Let C be an m x p Cauchon diagram. We say that a sequence 

((io,ao), (^i,ai), («*,«*)) 
is a lacunary sequence with respect to C if the following conditions hold: 
1. t > 0; 
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2. the boxes (ii,ai), (i2,«2), (^t,at) are white in C; 

3. 1 < if) < ii < ■ ■ ■ < it < m and 1 < ao < ai < ■ ■ ■ < < p; 

4. If it < i < m and at < a < p, then (z, a) is a black box in C; 

5. Let s G {0, 1}. Then: 

• either (z, a) is a black box in C for all ig < i < is+i and < a, 

• or {i, a) is a black box in C for all ig < i < is+i and ao < a < a^+i; 

(See Figure [2] where := Xs+i and Xk := {ik,ak)-) 

6. Let s G {0, 1}. Then: 

• either {i, a) is a black box in C for all is <i and < a < a^+i, 

• or (z, a) is a black box in C for all i < ig+i and < a < a^+i. 

(See Figure [31) 
Example 3.2. Consider the 3x3 Cauchon diagram 







m 


□ 






u 







It is easy to check that the sequence ((1, 2), (3, 3)) is a lacunary sequence (with respect to 
C). 

In [2], lacunary sequences were studied in order to prove the normal separation of the 
prime spectrum of the algebra of m x p quantum matrices. Moreover, Cauchon studied 
the existence of such sequences. In particular, he proved the following result. 

Lemma 3.3. Let (j, (3) G |1, m] x Assume there exists jo > j and f3o > (3 such that 

(jO)/3o) is a white box in C. Then there exists (/,7) G x such that 

1. I > i and 7 > (3; 

2. the box (/, 7) is white in C ; 

3. • Either {i, a) is a black box in C for all j < i < I and (3 < a, 
• or {i, a) is a black box in C for all j < i < I and a < 7; 
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Figure 2: Condition 5 of Definition 13.11 
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Figure 3: Condition 6 of Definition 13.11 



4- • Either {i, a) is a black box in C for all j < i and (3 < a < ^, 
• or {i, a) is a black box in C for all i <l and /3 < a < 7. 

Proof. We give a proof of tliis lemma for the convenience of the reader. We distinguish 
between three cases. 

Case 1: The first /3 columns of the Cauchon diagram C' obtained from C by deleting 
the first j rows are black. By assumption, there is a white box in . Let 7 be the smallest 
column- index of a non-black column in C*. Hence 7 > /3. We know that there is a white 
box in the column 7 of C^, so 

/ := min{/c \ j <k and the box (^,7) is white} 

exists and j < I < m. 



9 



We claim that (^,7) satisfies the properties of the lemma. Indeed, if j < i < I and 
q; < 7, then the box {i, a) is black as by construction the column indexed a of is black. 
So the second assertion of 3. is satisfied. 

Similarly, if j < i and ^ < a < 7, then the box {i, a) is black, so that the first assertion 
of 4. is satisfied. 

Case 2: The first j rows of the Cauchon diagram obtained from C by deleting the 
first /3 columns are black. By assumption, there is a white box in C^. Let / be the smallest 
row-index of a non-black row in C^. Hence I > j. We know that there is a white box in 
the row I of C"^, so 

7 := min{r7 \ (5 < r] and the box (/, rj) is white} 

exists and (3 < •y < p. 

We claim that (/, 7) satisfies the properties of the lemma. Indeed, if j < i < / and 
/3 < a, then the box {i, a) is black as by construction the row indexed i of is black. So 
the first assertion of 3. is satisfied. 

Similarly, if i < / and (3 < a < j, then the box (i, a) is black, so that the second 
assertion of 4. is satisfied. 

Case 3: The first /3 columns of the Cauchon diagram obtained from C by deleting 
the first j rows are not entirely black, and the first j rows of the Cauchon diagram 
obtained from C by deleting the first (3 columns are not entirely black. 

We set 

k :— min{s \ j < s and there exists 1 < 77 < /3 such that the box (s, 77) is white} 

and 

7 := mm{r] \ (3 < rj and there exists 1 < s < j such that the box (s, rj) is white}. 

We have j < k < m and (3 < j < p. 

Observe that the box {k, 7) is also white because C is a Cauchon diagram. Now we can 

set 

/ := min{s \ j < s and the box (s, 7) is white}. 

We claim that (/, 7) satisfies the properties of the lemma. Indeed, if j < i < I and 
q; < 7, then the box (i, a) is black. Otherwise, {i, a) is white and by construction of 7 
there exists 1 < s < j such that the box (s, 7) is white. As C is a Cauchon diagram, this 
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In this picture, the area in red cor- 
responds to the Cauchon diagram 
C", and the box marked with a W 
is white. 



Figure 4: Proof of Lemma [3.3[ Case 1. 
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In this picture, the area in red cor- 
responds to the Cauchon diagram 
C^, and the box marked with a W 
is white. 



Figure 5: Proof of Lemma [3.3[ Case 2. 




In this picture, the areas in red cor- 
respond to the two containing 
a white box, and the boxes marked 
with W are white. 



Figure 6: Proof of Lemma [3.3[ Case 3. 
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would force {i, 7) to be white, contradicting the minimahty of I. Hence, {i, a) is black, and 
the second assertion of 3. is satisfied. 

Now, assume that i < I and ^ < a < 7. We distinguish between two cases to show 
that the box (i, a) is black, li j < i, then the same argument as above shows that (i,Q;) 
is black. Indeed, otherwise {i,a) is white and by construction of 7 there exists 1 < s < j 
such that the box (s, 7) is white. As C is a Cauchon diagram, this would force (i, 7) to be 
white, contradicting the minimality of Hence, {i, a) is black. Finally, if i < j, then by 
construction of 7, the box (i, a) is black. Hence, the second assertion of 4. is satisfied. □ 

Fix (j, /9) G |l,m] X The proof of the previous lemma actually gives us an 

algorithmic way to produce a lacunary sequence starting at {j, . 

Algorithm 1 Algorithm constructing a lacunary sequence with respect to C from any box 
input: Let C be a Cauchon diagram, and fix (j, /3) G x 

output: a lacunary sequence ((io, ao), (^1, Q!i), (^t, ctt)) with respect to C such that 

(io,ao) = 

start: io := j, ctQ := (3, t := 0, i := it, a := at 

while the Cauchon diagram Cj^^ deduced from C by deleting the first i rows and the 
first a columns is not completely black do 

if the first a columns of the Cauchon diagram C obtained from C by deleting the 
first i rows are black then 

at+i :=the smallest column-index of a non-black column in 
and it+i :— min{/c \ i < k and the box {k, at+i) is white} 

else 

if the first i rows of the Cauchon diagram obtained from C by deleting the 
first a columns are black then 

if+i :— the smallest row-index of a non-black row in C"' 
and at+i :— m.m{r) | a < 77 and the box {it+i, rj) is white} 

else 

ctt+i := min{77 1 a < 77 and there exists 1 < s < i such that the box (s, 77) is white} 
and it+i :— min{s \ i < s and the box (s, at+i) is white} 
end if 
end if 
t:=t + l. 
end while 

return ((io, ag), (ii, cti), {it, at)). 
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Corollary 3.4. Let {j,f3) G [l,""!] x Then there exists a lacunary sequence 

((io,ao), («i,ai), i}t,Cit)) 

with respect to C with 

(io,ao) = 

Proof. Algorithm [1] produces such a lacunary sequence thanks to the previous lemma. □ 
Example 3.5. Consider the Cauchon diagram 



Applying the above algorithm starting from the box (1, 2), we obtain the lacunary sequence 
((1,2), (2,3)). Note that this is not the only lacunary sequence starting at (1,2). Indeed, 
one can easily check that ((1, 2), (3, 3)) is another lacunary sequence starting at (1,2). 

Remark 3.6. If j = m or /3 = p, then the only lacunary sequence (with respect to any 
m X p Cauchon diagram) is ((j, /3)). 

4 Recognition 

Fix a Cauchon diagram C. If {j,f3) G |l,'ml x |l,p], then we denote by LCj^/j the set 
of all lacunary sequences (for C) that starts at (j, /3). The previous section shows that 
LCj^/s ^ 0. If (j,^) := ((jo,/3o), (ji,/3i), {jt, f^t)) G LC^-^ is a lacunary sequence starting 
at {j,/3), then we set: 

Aj^ := [jo, ji, ...,jt \ l3o,(3i,...,f3t]eO (MmAC)) ■ 

Proposition 4.1. Let M G Aim,pO^)- For all (j, /?) G [l,""^! x |l,p], fix a lacunary se- 
quence (j,^) = {{jo, (3o), {ji, f3i), {jd, (3d)) G LCj^i3, and set Af^^ := A^^. Assume that 
either (a) or (b) below is true. 

(a) The matrix M is tnn and belongs to the tnn cell Sq associated to C. 

(b) For all (j, P), we have Af^M) = if (j, /3) G C and Af^^M) > z/ (j, /3) ^ C. 
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Then, 

for all lacunary sequences (i, a) = ((^o; c^o)) (^i; cti); {hi 

Proof. We prove by a decreasing induction on (j, /3) G fl,"^] x (with respect to the 

lexicographic order) that 

(i) For all (/, 5) > (j, then ti^s > and U^s = if and only if (/, 5) G C; 

(ii) For all lacunary sequences (i, a) = {{iQ,ao), {ii,ai), {it,at)) G LC-,-,^, and for all 
{k,ri) G -E" with {k,ri) > {j,f3), we have 

where < {k,r]) < {is+i,as+i). 

(By convention, if s = — 1 (ie if {k, rf) = (j, /3) = (zq, Oo)), then we set 

[to,...,ts I ao,...,«.](M('='^)) = l. 

On the other hand, if s = t, then we set {it+i, at+i) = {m + l,p + 1).) 

Note that condition (i) is automatically satisfied under assumption (a) thanks to The- 
orem 12.71 

As observed in Remark l3.6[ the only lacunary sequence starting at {m, (3) is ((m, /?)). 
So it is clear that A^^(M) = A^^^(M) = Xm,i3 = tm^p for all lacunary sequences starting 
at (m, This proves (ii) in the case where j = m. Moreover, under assumption (b), 
as A^^^(M) = if (j,/3) G C and A5;^(M) > if (j, /3) ^ C, we obtain that tm,i3 = if 
(m, /3) G C and t^,/? > if (m, /3) ^ C, which proves (i) in the case where j = m. 

Let j < m, and assume (i) and (ii) are true for all (j',/3') > (j, /3). In particular, we 
assume that ti^s > for all (Z, 6) > (j, (3) with equality exactly when (/, 6) G C. 

To prove (ii) at step {j,f3), we proceed by a decreasing secondary induction on {k,ri). 

When {k,ri) = (m + l,p), the result is trivial. So let (j, /3) < {k,ri) < (m + l,p) 
and assume the result established when {k,ri) is replaced by {k,ri)~^. In particular, for the 
lacunary sequence (i, a) = {{io,ao), (ii,ai), {it, at)) G we assume that 

^ia{M) = [^o, ...,^s\ao,.. ■ , a,](M('=''')^)ti,^,,„,^, ■ ■ ■ti,,^,, 

where (z^,^^) < {k,r])+ < {is+i,as+i). 
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For each {k,ri), exactly one of the following five cases must occur 
Case 1: {k,r]) = (is,"^); 
Case 2: k = ig and rj > Og] 

Case 3: is < k < is+i and r/ G {1, . . . , ao} U {ai, . . . , as}; 
Case 4: ig < k < z^+i and r] > a^+i; 

Case 5: is < k < is+i and there exists h G [0, s] such that ah < rj < ah+i- 
We investigate each case separately to prove (ii). 

• Case 1: {k,r]) = {is,as). 
If s = 0, then we have 

[zo, ...,is\ao,.. . ,as]iM^'''''^^) = U^^ao = 1 ■ Uo,ao, 

and the result is proved. 

Assume now that s > 0, so the box {is,as) is white in C. As = {k,r]) > 

(j, (3) = {io, cto), we have ti^^^s > by the primary induction hypothesis on (j, (3). From [U 
Proposition 3.7], we deduce: 

[^0, . . . , I ao, • • • , as]{M^'''''^^) = [to, | «o, • • • , «.-i](M(^'°^)) ■ U^^^^, 

as required. 

• Case 2: k = ig and rj > Og- In this case, it follows from [8| Proposition 3.11] that 

[io, ...,tg\ao,..., a,](M(^''')+) = [zq, . . . , | ao, ■ • ■ , a.](M('=''')), 

as desired. 

• Case 3: ig < k < ig+i and 77 G {1, ... , ao} U {ai, . . . , a^}. In this case, it follows from 
[H Proposition 3.11] that 

[to,...,ts I ao,...,a.](M(^'')^) = [^o,...,^. | ao, • • • , «.](M('=''')), 

as desired. 

• Case 4: ig < k < ig+i and r/ > a^+i. 

If (A;, 77) is black in C, then tk,rj = by the primary induction hypothesis on (j, (3), and 
so it follows from [HI Proposition 3.11] that 

[to,...,ts I «o,...,as](M(^'''') = [^o,...,^. I «o,...,«.](M(^'^)^), 
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as desired. 

Now we can assume that {k, rf) is white. In this case, it follows from [HI Proposition 
3.13] that the difference 

[io,...,z. I ao,...,a.](M('='''^^)-[zo,...,z. | ao, • • • , a.](M('='''^) 

is a scalar multiple of tk^a^- However, as (/c, rf) is white and the sequence ((^o, ao); (^i? Q^i); {.Hi 
is lacunary, we deduce from the fifth condition of the definition of lacunary sequences that 
the boxes (/c, a) are all black for < a < a^+i. In particular, (fc, a^) is black in C, and 
so tk,as = by the primary induction hypothesis on (j, (3). Hence, we get 

[io, . . . , I ao, . . . , ttsllM^'^'^)^) = [zo, . . . , I ao, . . . , ttsJlM^'^'")), 

as required. 

• Case 5: ig < k < i^+i and there exists h G |0, s] such that ah < r] < ah+i- If {k,ri) 
is black, then tfc ,, = by the induction hypothesis on (j, and so it follows from [HJ 
Proposition 3.11] that 

[io,...,is I «o,...,«s](M('^'^)) = [to,...,ts I «o,...,«.](M(^'^)^), 

as desired. 

From now on, we assume that {k, rf) is white, so that t^^^ > (by the primary induction 
hypothesis on (j, 

First, as {k^rf) is white with k > ih and at < f] < a^+i, the first part of the sixth 
condition of a lacunary sequence is not satisfied for the lacunary sequence (i, a). Hence, 
the second part has to be satisfied, that is: {i, 6) is black for all i < ih+i and ah < S < ah+i- 

We claim that the sequence {{i^rfj, {ih+i,ah+i), ■ ■ ■ , {it,at)) is lacunary for every i < 
ih+i- For, we just need to prove that {i',a') is black when i' < ih+i and rj < a' < ah+i 
as this will prove that the second parts of the fifth and sixth conditions are satisfied for 
the first step of the lacunary sequence. Under these conditions, we have i' < ih+i and 
O-h < V c^' < o.h+ii so the above shows that the box {i\ a') is indeed black. In particular, 
note for later use that (i, rfj is black for all i < ih+i- 

Let i G {io, . . . ,ih}- Note that {i,ri) > (io^^o) = {j,f3)', so (i) and (ii) apply to any 
lacunary sequence starting at {i,ri), by the primary inductive hypothesis. We apply (ii) to 
the lacunary sequence ((i, rj), {ih+i, • • • , {it, c^t)) in two instances. Recall that h < s; 
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so {is+i,as+i) is in this lacunary sequence. Moreover, as i < is < k < is+i, we have 
{is,as) < {k,ah) < {is+i,(ys+i) if h < s and {i,r]) < {k,ah) < {is+i,as+i) if h = s; so, in 
both cases, we get 

[i,ih+i, ■■■,it I V,Oih+i, ■ ■ . ,at](M) = 

[i,ih+u...,is I ?7,ah+i,...,as](M(''''"'^)) -tj.+i.a.+i • ... ■ 

Next, apply (ii) at the pivot {i, rf) (so we are in the case s = —1) to obtain 
[i,ih+i, I V,Oih+i, ■ ■ .,at]{M) = l.t^^^ ■ U^^^^c^h+i ' ■■■ ' 

Thus, 

Recall that, for all i G {iq, . . . , ih}, the box (z, rj) is black. Hence, tj^^ = by the primary 
induction hypothesis on (j, /3) (which we can apply as {i.rj) > {io,v) > (^O)Cto) = {hP))- 
So we get 

[i, . . . , r], ah+i, . . . , as](M('''"'^)) ■ ■ ... ■ U^^o,, = 0. 

Similarly, as the boxes {is+i,as+i), iit^oit) are all white, we get that ti^^^^aa+n ■ ■ ■ i'tit,at 
are all positive, so 

I r/, «;,+!,..., a,] (M('='°'')) =0, (1) 

for all i G {iq, . . . , z/i}. 

Note now that the sequence ((«/i+i, Oih+i), • • • , (^t, is lacunary. So following the same 
reasoning as above we get that 

[4+1, ...,is\ ah+i, . . . , a,](M('^'°'^)) = • • • K.a, > 0. (2) 

We then deduce from ([T]), ([2]) and Sylvester's identity that 

[io, • • • 4, 4+1, • • • , 4 I "0, • • • , tth-i, ^, tth+i, • • • , as]iM^'''°"'^) = 0. 

Hence, we deduce from |8i, Proposition 3.13] that 

[zo, ...,Zs\ao,..., a,](M('=''')) = [^o, . . . , | ao, • • • , «.](M('=''')^), 

as desired. 
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Hence, we have just proved that (ii) is true at step (j, It remains to estabhsh (i); 
that is we need to prove that tj^^ > with equahty if and only if (j, f3) G C. Note that 
this is automatic under assumption (a) thanks to Theorem 12. 71 Under assumption (b), by 
(ii) at step (j, /3) with (/c, r/) = (j, we have 

As the boxes [ji, Pi) with 1 < I < d are white, we deduce from the induction hypothesis that 
> f*^^ ^- '-^fA^^ — ^ "with equality exactly when (j, /3) G C (by hypothesis), 
we get that tj^^ > with equality exactly when (j, /3) G C This proves (i). □ 

The following two results now follow easily. 



Proposition 4.2. Let M G Airn,p{^)- We keep the notation of Proposition Assume 



that M is tnn and belongs to the tnn cell associated to C. Let (j, /3) G [l,""!] x 
Then we have ^lp{M) = if (j, /3) G C and A5;^(M) > if (j, P)(^C. 

Proof. Note that A^^(M) > as M is tnn. This result follows easily from the previous 

5^ associaieu Lu o, we iiavt; i^j^^f}^ 
s > by Theorem 12. 7[ Hence, A^^(M) = if and only if tj^^ = tj^^^jj^ = 0; that is, if and 
only if (j, /3) G C as M belongs to the tnn cell associated to C. □ 



proposition. Indeed, A^^(M) = tjo,/3o ' ' ' ''^jd^Pd- boxes {js,(3s) for s > are 

white in C and M belongs to the tnn cell 5*^ associated to C, we have t,-^ « > for all 



Proposition 4.3. Let M G AimA^)- keep the notation of Proposition Assume 
that, for all (j, /3), we have ^lp{M) = if (j, /3) G C and Af^p{M) > if (j, 13) ^C. Then 
M is tnn and belongs to the tnn cell associated to C . 

Proof. By Proposition 14. we know that each ti^a > and that ti^a > precisely when 
{i,a) G C. The result now follows from Theorem 12.71 □ 

Propositions 14.21 and 14.31 together provide a test for membership of a fixed totally 
nonnegative cell. Each cell is specified by a Cauchon diagram. For a given m x p Cauchon 
diagram C assume that a lacunary sequence has been chosen for each box (j, /3) and for a 
given real mxp matrix M check whether A^^(M) > for each white box and A^^(M) = 
for each black box. Then M is totally nonnegative and is in the cell specifed by C if and 
only if these mp statements are all true. This is summarised in our main result. 

Theorem 4.4. Let C be an m x p Cauchon diagram. For all {j,f3) G 11,^.] x |1,J9], fix 
a lacunary sequence (j,^) = ((jo,/3o), (ji, /^i), {jd, (3d)) e LCj-/3, and set Aj;^ := A^^. 
Let M G AimA^)- following are equivalent. 
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1. M is tnn and belongs to the tnn cell associated to C . 



2. For all (j, /3) G [l,m] x we have Af[p{M) = if {j,f3) G C and Af[p{M) > 

Note that this test only involves mp minors. In the case that C is the Cauchon diagram 
with all boxes coloured white, the test states that a real matrix M is totally positive if and 
only if each final minor of M is strictly positive. As discussed earlier, this is the well-known 
Gasca and Pena test, but applied to final minors rather than initial minors. 

Example 4.5. Consider the following Cauchon diagram: 



C 







■ 




□ 











Then a 3 x 3 real matrix M is tnn and belongs to the cell associated to C if and only if 
Ai,i(M) = [13|12](M) > Ai,2(M) = [12|23](M) > Ai,3(M) = [1|3](M) = 

A2,i(M) = [23|12](M) = A2,2(M) = [23|23](M) = A2,3(M) = [2|3](M) > 

A3,i(M) = [3|1](M) > A3,2(M) = [3|2](M) > A3,3(M) = [3|3](M) > 0. 

(The lacunary sequences have all been obtained by using Algorithm [H) 

/ 16 5 \ 



It is easy to check that the matrix M 



12 6 3 



satisfies the above nine condi- 



V 4 2 ly 

tions. Hence, we deduce from Theorem 14.41 that M is tnn and belongs to the tnn cell Sq 
associated to C. 
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